The Universality Class of Diffusion Limited Aggregation and Viscous Fingering 
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We investigate whether fractal viscous fingering and difltusion hmited aggregates are in the same 
scahng universality class. We bring together the largest available viscous fingering patterns and a 
novel technique for obtaining the conformal map from the unit circle to an arbitrary singly connected 
domain in two dimensions. These two Laplacian fractals appear different to the eye; in addition, 
viscous fingering is grown in parallel and the aggregates by a serial algorithm. Nevertheless, the 
data strongly indicate that these two fractal growth patterns are in the same universality class. 



Laplacian fractals are paradigmatic examples for the 
spontaneous growth of fractal patterns in natural sys- 
tems. In particular, two such examples have attracted 
an enormous amount of interest: viscous fingering and 
diffusion limited aggregation (DLA). 

Viscous fingering is realized when a less viscous fluid 
displaces a more viscous fluid contained in the narrow 
gap between two glass plates (a Hele-Shaw cell). When 
the less viscous fluid is inserted through an aperture in 
one of the glass plates, a pattern is formed, as shown in 
Fig. n .2] . The fluid velocity v of the displaced fluid in a 
Hele-Shaw cell satisfies Darcy's law, v cx Vp, where p is 
the pressure. To a very good approximation, the viscous 
fluid is incompressible (i.e., V • v = 0), so the pressure in 
the viscous fluid satisfies the Laplace equation V^p = 0. 

DLA is realized 3] as a computer experiment in which 
a fractal cluster is grown by releasing a fixed size random 
walker from infinity, allowing it to walk until it hits any 
particle already belonging to the cluster. Since the parti- 
cles are released one by one and may take a long time to 
hit the cluster, the probability field is stationary and, in 
the complement of the cluster, we have again V^p — 0. 





FIG. 1: Digitized image of an experimental viscous fingering 
pattern. Air (black) is injected into the oil-filled gap (white). 
The pattern is approximately 22 cm in diameter, cf. 



FIG. 2: A diffusion hmited aggregate with 50,000,000 parti- 
cles. Image courtesy of EUak Somfai. 



A typical DLA cluster is shown in Fig. |21 

Mathematically, these problems are similar but not 
identical. In both cases, one solves the Laplace equa- 
tion with the same boundary condition at infinity, i.e. 
Vp = const X r/r as r oo. However, in viscous fin- 
gering, each point on the fractal's boundary advances at 
a rate proportional to Vp, whereas the DLA accretes 
one particle at a time, changing the Laplacian field after 
each such growth step. Thus, we refer to viscous finger- 
ing and DLA as parallel and serial processes, respectively. 
In addition, the ultraviolet regularization differs; in DLA, 
p ~ on the cluster and the regularization is provided 
by the particle size. In viscous fingering, one solves the 
problem with the boundary condition p — an where a is 
the surface tension and k is the local curvature. Finally, 
viscous fingers are grown in a finite gap and are not truly 
2-dimensional. Accordingly, one can ask whether these 
two fractal growth problem are in the same scaling uni- 
versality class. But to answer this question, one must 
first define what one means by a "scaling universality 
class." 

Definition of the scaling universality class: The 

first test of correspondence between the fractal patterns 
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has to do with their fractal dimension Dq. Denote by 
Rn the radius of the minimal circle that contains a frac- 
tal pattern. The fractal dimension is defined by how the 
mass Mn contained within this circle (number of parti- 
cles for DLA, area for viscous fingers) scales with i?„, 
Mn ~ Rn" ■ Measurements of this type indicated a value 
Dq « 1.71 for both problems, motivating many authors 
to express the opinion that these two problems are in the 
same universality class 0, ■ Obviously, the fractal di- 
mension by itself is not sufficient, and a more stringent 
definition is necessary. 

We propose here that the identity of the scaling prop- 
erties of the harmonic measure is a sufficient test for two 
Laplacian growth problems to be in the same universality 
class. The harmonic measure is the probability for a ran- 
dom walker to hit the boundary of the fractal pattern. It 
determines the growth in both problems, being propor- 
tional to Vp at the boundary. Suppose that we know the 
probability measure /i(s)ds for a random walker to hit an 
infinitesimal arclength on the fractal boundary. We com- 
pute the probability Pi{e) = J^'^j^ box '^('^)^'*' then 
define the generalized dimensions 5] via 



logV^^f^ PHe) 
D„ = lim ' ^ ' 

e^o [q - 1) loge 



(1) 



where N{e) is the number of boxes of size e that cover 
the fractal boundary. The g > branch of this func- 
tion probes the high probability region of the measure, 
whereas the q < Q branch stresses the low probabilities. 
For q — > O"*" we find the fractal dimension Dq. An equiv- 
alently useful description j3| is provided by the scaling 
indices a that determine how the measure becomes sin- 
gular, Pi{e) ^ e", as e ^ 0. This set of indices is ac- 
companied by /(a) which is the fractal dimension of the 
subset of singularities of scaling exponent a. The rela- 
tion of these objects to Dq is in the form of the Legendre 
transform 



a{q) 



d[{q-l)Dq] 



f{a) =qa{q) - {q - 1) Dq 



(2) 
(3) 



In fractal measures defined on simple fractals, the values 
of a are usually bounded from above and from below, 
Q^mm < < ctmax- When the fractal measure fails to 
exhibit power law scaling everywhere (say, Pi{e) ^ e^*^ 
somewhere), one finds a phase transition in this formal- 
ism with one of the edge values Umin or Umax ceasing 
to exist. When the spectrum of exponents a and their 
frequency of occurrence f{a) in two problems are the 
same, we refer to these problems as being in the same 
universality class. We note that a similar criterion was 
used to state that different dynamical systems are in the 
same universality class at their transition to chaos 
but not many other physical problems yielded to such a 
stringent test, simply because it is not easy to compute 
with enough precision the scaling properties of fractal 
measures. For DLA, such an accurate computation was 




FIG. 3: Conformal map of the unit circle aimed to fit the pat- 
tern in Fig. Q The map was constructed using the algorithm 
proposed in the text, transforming the unit disc represented 
by the black area. This figure demonstrates the accuracy of 
our method for reproducing complex patterns. 



achieved 0; in this Letter, we report also an accurate 
computation for experimental viscous fingering patterns, 
allowing us to test whether the two problems are in the 
same universality class. 

The apparent complexity of the viscous fingers is dis- 
couraging for any attempt to calculate the harmonic mea- 
sure reliably. The naive method for computing the har- 
monic measure for DLA is to "probe" the interface with 
random walkers and perform frequency statistics 
While the statistics at the outer tips is reasonable 
the deep fjords are visited extremely rarely by the ran- 
dom walkers. Similarly, in viscous fingering experiments, 
the measure is usually estimated from the velocity of the 
interface ^2 ■ Obviously, this reveals the outer tips which 
move with an appreciable velocity, but leaves unmeasured 
the harmonic measure on the fjords, which almost do not 
move at all. One has to look for alternative methods. 

The harmonic measure from iterated conformal 
maps: The Rieniann theorem guarantees that there ex- 
ists a conformal map from the exterior of the unit circle 
u! = e*^ to the exterior of our viscous fingering patterns. 
Having such a conformal map, say ^{uj), the harmonic 
measure is simply obtained as 1/|$'((jj)| up to normaliza- 
tion. The actual construction of such a map for a given 
well-developed fractal pattern is, however, far from obvi- 
ous, and it has never been accomplished before. We will 
demonstrate now that the method of iterated conformal 
maps 13, Ijlj can be adapted to solve this problem in a 
very efficient way. We use data acquired in experiments, 
in which air is injected into oil from a central orifice. The 
oil is confined between two circular glass plates, 288 mm 
in diameter and each 60 mm thick 0. The plates are 
separated by 0.127 mm and are flat to 0.13 /.tm (optically 
polished to one-quarter of a wavelength). An oil buffer 
surrounds the plates. The data presented here were ob- 
tained with silicone oil with a dynamic viscosity fi = 345 
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mPa-s and surface tension cr=20.9 mN/m at 24°C; addi- 
tional experiments were conducted with silicone oil with 
a dynamic viscosity fi = 50.8 mPa-s and surface tension 
(7=20.6 mN/m at 24°C. The pressure difference Ap be- 
tween the injected air and the oil buffer ranged from 0.1 
to 1.25 atm. The 1 mm diameter central hole through 
which the air is injected cannot be seen in Fig. ^ rather 
the 22 mm diameter central solid circular region (black) 
is a mask introduced after the images have been digi- 
tized. The interface near the injection point could not be 
extracted because of light contrast limitations. 

The images are digitized and each pixel is given an in- 
dex j from 1 to A''. As a first step, we rescale all the pix- 
els and apply a uniform shift such that the central black 
region becomes the unit circle. Next, we record the po- 
sition of each pixel in the complex plane, denoting it by 
Zj = Xj +iyj. Next, we sort the pixels in ascending order 
with respect to \zj\, which is the minimal distance Vj of 
the j'th pixel from the interface of the fractal pattern. 
We divide the sorted pixels into m groups, such that the 
fc'th group contains the pixels which have absolute values 
smaller than those in the (/c-f l)'th group. Within the in- 
dividual groups, we once again divide the pixels into two 
subgroups, where the first subgroup contains the pixels 
with associated radii rj larger than a threshold value 6. 
Typically, we choose S to be five times the pixel size. We 
build the conformal mapping from the unit circle to the 
viscous fingers by expanding stepwise the unit circle in 
the interior part of the viscous fingers. The size of the ex- 
pansion is given by rj . In other words, in the n'th step we 
build a conformal mapping (w) from the uj complex 
plane to the z complex plane. At each step, we employ an 
auxiliary mapping (pg^\ that locally perturbs the unit cir- 
cle with a semicircular bump centered at the position e*^ 
and of linear size a/A. The construction is done such that 
we first expand around the points in the first subgroup of 
the first group, then we proceed to the second subgroup. 
Once we have finished expanding around the first group, 
we move on to the second group, etc. Most pixels in the 
vicinity of the pixel Zj will be covered when expanding 
the conformal mapping around Zj] therefore, they can be 
disregarded in all future iterations. In fact, a majority of 
the pixels will be covered by expansions around nearby 
points. Thus we need to perform a significantly smaller 
number of expansions than the number of pixels. The 
aforementioned division of the pixels is a simple way to 
speed up our method, i.e. to reduce the number of nec- 
essary expansions. We first expand the mapping around 
the points x'^^Hzj) with large radii (rj > 6) and close to 
the unit circle. The conformal mapping is constructed 
by successive iterations. In each iteration step, the map- 
ping covers more and more of the viscous fingers, and the 
number of pixels not being covered decreases. 

Assume that we have reached the n'th iteration step. 
Considering the pixels that are not covered, we now use 

the inverse mapping, x'"' = , to map these pix- 

els Zj to the exterior of the unit circle, i.e. we map them 
from the z-plane to the w-plane. The size -y/ A„+i of the 
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FIG. 4: Log-log plots of ^ii^) the box size e for 

selected values of q. We find reliable power law scaling for 
q > -0.25. 
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FIG. 5: /(a) curves for the harmonic measure of the exper- 
imental viscous fingers patterns and its comparison with the 
/(a) curve of the harmonic measure of DLA. (a) is averaged 
over 15 patterns obtained with a silicone oil with a dynamic 
viscosity fi — 345mPa-s and (b) is averaged over 15 patterns 
with p — 50.8mPa-s. 



(n + l)'th bump to be centered at the position Zj is taken 
to be 



$'(n)(eie„+i) 



■, where On+i ^ aigx^^H^])- 

(4) 

From the unit circle, we now advance along the shape by 
expanding around the next Zj, according to the scheme 
described above. 

The efficiency of this method is surprising; the number 
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FIG. 6: Close-up of the /(a) curves presented in Fig. |^ to- 
gether with the same curve computed from a samphng with 
random walkers. 

of iterations needed to build the conformal map remains 
rather low even for a complex pattern like that in Fig.^ 
Figure 121 presents the result of our method using the dig- 
italized image of viscous fingers shown in Fig.^ Indeed, 
the conformal map accurately reproduces the image. In 
general, the number of iterations is around 10000-20000 
for the maps that we analyze in this Letter. 

Results: With the conformal mapping at hand, we 
have the harmonic measure as One has to 

be cautious in evaluating the derivative, since inside the 
fjords, the measure becomes so small that it is difficult 
to resolve the fjords on the unit circle. The way to eval- 
uate the derivative involves a careful book keeping of the 
positions of the individual bumps added in each iteration 
step 01 . In Fig. ^ we present log-log plots of the sum 



EiLiP'i^) vs. the box size e for selected values of q. 
The estimated slopes, r(g) = {q— ^)Dq, are very reliable 
for positive values oi q. As seen in Fig. 21 for values of 
q < —0.25, the power-law scaling is lost. In DLA, the 
deep fjords become dominant for q < —0.2, and lead to a 
nonanalyticity in the f{a) curve, a phenomenon known 
as a "phase transition" in the thermodynamic formalism 
of the fractal measure 0| ■ 

We can now use the measured values of T(q) to ex- 
tract the f(a) curve. The result is shown in Fig. [SI 
together with the same curve for the harmonic measure 
of DLA. While there is no alternative way to compute 
this function in its entirety, we can test the accuracy of 
our results by looking at /(a) for g > 1, which charac- 
terizes the highly probable regions of harmonic measure 
near the tips. We sampled the tips of the experimental 
data with random walkers, a good method as long as one 
studies the part of the multifractal spectrum for q > 1. 
For the sampling, we use 10^ random walkers. Increasing 
this number did not lead to a significant change in the 
results. Figure ini shows a close-up of the /(a) curve pre- 
sented in Fig. [31 alongside the one computed from the 
random walk samples. There is almost no difference in 
the left hand branches, and we thus conclude that our 
method is fully adequate for representing the tips of the 
viscous fingers. The estimate of amm is the same (within 
the uncertainty of ±0.06) as the one computed for DLA, 
which was amm = 0.67 |l5l |. 

In conclusion, we have presented a powerful method for 
the construction of a conformal map from the unit circle 
to complex fractal patterns. Using this conformal map, 
we can compute accurately the fractal measure and de- 
termine its scaling properties. The comparisons shown in 
Figs. OandElstrongly indicate that viscous fingering and 
DLA are in the same universality class within the experi- 
mental uncertainty. Future study can determine whether 
this correspondence extends to the nonanalyticity in the 
/(a) curve. 
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